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This paper presents a solution to the nonlinear small x "projectile side" evolution equations as derived by Balitskii 
in 1996. The solution is based on functional Fokker-Planck methods. The fixed point at small x is explicitly 
calculated and all correlation functions in this limit are determined. They show clear saturation and unitarization 
properties. Scaling laws that hold during the saturation phase and throughout the whole course of the evolution 
are established. The corresponding Langevin equations are given as a basis for numerical simulations opening 
the field for future studies of dynamical issues of the evolution not analytically accessible. The methods used 
may be extended to the "target side" equations of Jalilian-Marian, Kovner, Leonidov and Weigert. 



1 Introduction 

Recent years have seen a growing activity centered around systems and experiments involving large numbers 
of partons. Experimental activities contributing to this newly focussed interest have been conducted and 
arc planned with widely varying focus such as the deep inelastic scattering (DIS) experiments at HERA, the 
heavy ion experiments at CERN with LHC and RHIC at BNL. One common denominator is that at all these 
experiments, be it because of high energies in ep collisions or because of intrinsically higher numbers of par- 
ticipants in eA and AA collisions, QCD effects are dominated by the large numbers of gluons involved. This 
fact will clearly gain importance as the energy or the size of the participants grows with newer experiments 
at RHIC and LHC. 

As a consequence, one has to face up to the presence of these many degrees of freedom if one wants to 
understand the qualitative changes in the behavior of these systems as compared to lower energy and smaller 
targets. The main new feature expected is a "saturation of gluon densities" with implications for experiment 
ranging from global features such as unitarization effects in cross sections to their more subtle and harder to 
measure aspects in exclusive measurements. 

There have been many attempts to calculate, or at least to formulate the tools to do so, the energy dependence 
-in the guise of Bjorken x dependence- of cross sections and other observables as well as their A dependence. 
Characteristic for all these attempts is the basic fact that they, by necessity, use nonlinear features of QCD 
to incorporate high density effects. This activity has reached a point at which eA options for the HERA and 
RHIC programs are being actively considered. 

Still, besides many qualitative steps forward that all seem to confirm the generic ideas of "gluon saturation," 
the theoretical understanding remains limited. Partly this is due to the limited information presently available 
from experiment. The bulk of data used today is concerned with the total DIS cross sections and was measured 
at HERA. Although these are perfectly compatible with very tightly constrained saturation fits ||l|,|^, they 
are not fully conclusive, since conventional interpretations of the same data using Dokshitzcr-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution appear to be equally successful in reproducing the general trends in the 
data. This latter evolution is intrinsically linear and thus has no knowledge of the density effects important 
for unitarization and saturation. However, in comparing with experiment, the DGLAP machinery has been 
applied down to very small momentum transfers and serious conceptual criticism on the use of these 
equations so far beyond their region of validity is in order. Nevertheless, the fact remains that it has not 
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been conclusively shown that DGLAP evolution is violated in the present data. What one needs then is 
further experiments to provide the means to see DGLAP evolution fail and to conclusively demonstrate the 
presence of nonlinear effects, rather than just their compatibility with available data. The main effort on the 
theoretical side, however, has to focus on more quantitative treatments of nonlinear effects at small x and/or 
large A. 

Theoretical approaches which may form the basis for such endeavors have been formulated using widely 
different techniques. This often obscures their basic similarities and hinders comparison. To list but a few, 
nonlinear small x evolution equations have been presented by Mueller et al. , Balitskii |^,|^ , Kovchegov , 
Gribov et al. 1^, Levin et al. ||^-[lT|], Bartels and Ewerz as well as Jalilian-Marian, Kovner, Leonidov and 
Weigert (JKLW) |p^-|l7|. All these equations generalize the linear Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
equation for the gluon distribution to nonlinear -often infinite systems of coupled- equations for more general 
objects. They all show saturation and unitarization behavior at small x, although this often has only been 
demonstrated in simple limits such as the double logarithmic approximation. 

The relationships between many of these approaches has recently been investigated in which has finally 
put them on the same map. It had been known that the dipole approach of Mueller's as applied by Kovchegov 
to DIS was contained in Balitskii's approach to small x evolution. Ref. ||l8| has allowed to identify ||^,^ as 
a limiting case of the JKLW approach. This is not to say that the JKLW approach or Balitskii's equations 
have superseded the "simpler" evolution equations, quite to the contrary, there are clearly regimes in which 
the latter do capture the important physics. Being simpler, they will play their role in analyzing what is 
happening in new experiments. 

Moreover, although all these nonlinear equations have been received with interest, they have all mostly 
remained in the virgin state of a newly derived equation, without much knowledge about their solutions and 
concrete applications to physical problems. The more general and hence complicated the set of equations, 
the more this is true. The use of the more general approaches has up to now been -to some extent- more as 
a tool to organize simpler variants into genealogical trees than to provide new physics insight. It is simply 
the sheer complexity of the full nonlinear equations of Refs. ||^ and ||l^ [l^ that has hindered progress. To 
fully appreciate the scale of the problem, suffice it to say that they all result in infinite hierarchies of coupled 
equations for QCD correlators. Solving those is tantamount to solving a nonlinear field theory, albeit a bit 
simpler than QCD itself. However, the very way this "field theory" is presented has also precluded any 
numerical work that does not take recourse to truncations and other simplifications. Such approximations, 
however, would inherently limit the value of any results obtained, as it would at least partly eliminate the 
feature introduced with great care to capture the density effects, the nonlinearities themselves. 

This paper tries to tear down the brick wall which made progress come so slow. As it turns out, there exist 
methods to completely solve at least one of these infinite sets of evolution equations, the equations given in 
Ref. H]. As will become clear while doing so, the same methods should provide a good starting point to also 
tackle the somewhat more complex JKLW equations. 

In the early stages of rederiving the evolution equations of Ref. Q in preparation for jl^ it became clear 
that there is a much more compact way of presenting the results, quite analogous to the small x Wilson RG 
approach of JKLW, but with far less involved evolution kernels. It is this relative simplicity that makes these 
evolution equations the ideal candidate to search for further insight into their fixed point structure and to 
lay the groundwork for future numerical studies. The lessons learned in this example will undoubtedly prove 
useful in the further study of the JKLW equation, but will simultaneously give results for all the limiting 
cases already contained here. It will, in particular, provide a way to judge their respective regions of validity. 
The most prominent example in this case being dipole evolution as formulated by Ref. and applied to 
DIS in Ref. §. 

In analogy with the JKLW equation, Ref. j|]'s set of evolution equations can be written in form of a functional 
Fokkcr-Planck (FP) equation for a statistical weight that characterizes the expectation values of operators in 
the wave function of the target. This formulation provides the cornerstone for the progress reported here. It 
will lead to an understanding of the limiting behavior at small x and at the same time provide a formulation 
amenable to numerical treatment. 

To set the scene for these developments, Sec. || sketches a simple rederivation of Balitskii's infinite set 
of coupled RG equations for correlation functions at small x. In doing so a very brief discussion of the 
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relationships of the main nonhnear approaches so far in existence is provided. 

The rederivation itself has the purpose of emphasizing the basic physics ideas and concepts involved and 
provides a formulation in terms of a generating functional much more compact than the original version. 
While this in itself does not add anything to the physics results already known, it provides a convenient 
framework in which to extract all the known limiting or special cases and restricted versions of such evolution 
equations in a very transparent way. A short discussion on how to do this for a few select cases and of the 
physics implications of the approximations involved is given in Sec. |. 

Sec. ^ steps beyond the known results by converting the RG equation for the generating functional into an 
equation for a statistical weight, a formulation very similar to that used in the approach of JKLW [p^-[l7[ .f] 
Here emphasis is put on the fact that the resulting equation is of the form of a functional Fokker-Planck 
(FP) equation. This allows to interpret the statistical weight as a probability function(al) (in a very concrete 
sense) for the relevant operators at small x and to formulate the question of the existence and nature of a 
fixed point of the RG equation in terms of what one would normally call the equilibration properties of the 
FP operator involved. This indeed turns out to be a very powerful analogy that allows to identify the small 
X asymptotics of this evolution equation in full generality and at the same time paves the way for a new 
numerical approach. All subsequent results are based on this fact. It is obtained in two separate steps, the 
first, given in Sec. ^, is to find that there is no "restoring" force and hence one is dealing with a case of 
generalized Brownian motion. In a second step it is shown in Sec. ^ that this indeed corresponds to a case 
of ordinary Brownian motion for physical initial conditions allowing to eliminate some redundancy of the 
original presentation and to provide a yet simpler formulation. It is now a trivial consequence that there is 
a fixed point in this evolution and easy to show that it is attractive: the Brownian process asymptotically 
fills all available configuration space and one ends up with a very simple uniform statistical weight. In turn 
this means that in the approximation underlying this evolution, Ref . , QCD exhibits universality at small 
X. Cross sections and other measurable quantities become -in certain features- independent of the target 
at small enough x. This is a useful result even if it turns out that the approximations break down before 
the asymptotic limit is reached as it provides a clear overall tendency for the flow of observables within the 
region of applicability. More than that, the nature of the evolution kernel itself will allow to extract scaling 
laws that apply also far away from the fixed point and gives rise to predictions wherever there is a window 
of applicability of the evolution equations as such. 

Sec. is devoted to demonstrating one of the main new tools now available: The new approach allows to 
write down the Langevin equation corresponding to this FP/RG evolution. This on the one hand gives a very 
intuitive picture for the whole evolution -allowing to deduce scaling laws for the growth of the target and 
the cross section in the unitarization region. On the other hand one is then in a position to apply numerical 
simulations in order to study "dynamical" issues such as the rate of approach to the fixed point or, in other 
words, the practical aspects of how soon the system loses its memory of target specific scales and properties. 

Sec. 1^ discusses the asymptotic form of the correlators in the unitarization region, given some mild assump- 
tions about the negligibility of edge effects that should be good for sufficiently large targets. This exemplifies 
the scaling laws found in Sec. |^ and sheds some more light on the unitarization mechanism discussed in 
the previous sections. This also highlights some of the limitations in extracting quantitative information 
analytically that can be overcome by numerical simulations, such as the precise way evolution erases details 
of the initial conditions and reaches the asymptotic region. 

Sec. I will summarize the results and try to put them into context. Several appendices contain technical and 
notational odds and ends. 



2 A compact summary of Balitskii's equations 

The goal of this section is to set the stage for the fundamentally new results developed starting with Sec. |4| 
by summarizing the physics ideas behind the small x evolution mechanism underlying the results of jsj . The 
presentation in itself, however, is new and has not been given in this form before. It provides for a very 

^For a more thorough comparison of these two approaehcs see [H. 
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concise way of laying out the calculations involved and leads to the most compact formulation of the results 
within a single equation. I therefore hope this formulation will provide an easier access point to this field than 
other already published variants. It certainly makes it much easier to extract limiting cases and provides the 
simplest possible starting point for the studies laid out below. 

To begin with a physics motivation, I will take recourse to the parade ground example for saturation questions 
in small x physics, deep inelastic scattering (DIS) of virtual photons on hadronic targets of any size. The 
problem is set up in the frame in which the photon fluctuates into an energetic quark-antiquark pair long 
before it reaches the target, but where most of the energy resides in the target hadron which moves very fast. 
The scattering of the quark-antiquark pair is dominated by its interaction with the gluons in the target. Since 
the target hadron moves fast, the time evolution of the gluon fields is slowed by Lorcntz time dilation. Also, 
due to Lorcntz contraction, the gluon fields are well localized in the plane perpendicular to the direction of 
motion, which is taken to be the positive axis. The target can, therefore, be modeled by a distribution 
of static gluon fields localized at x^ = 0. As the scattering energy increases {x decreases) the gluon fields 
of the target change due to contributions of quantum fiuctuations. It is this evolution in x of the hadronic 
ensemble that one intends to describe in terms of the evolution equation. 

It proves convenient to use the light cone gauge = 0. In this gauge, loosely following Ref. Q, one takes 
the vector potentials representing the relevant gluon field configurations to be of the form 

A=b + dA with 0, b+ ^ P{x)S{xr) . (1) 

The question if this form is the only relevant one is in fact nontrivial and determines the region of validity of 
the present approach. This is discussed in detail in [T^ . 

The DIS structure function (which determines the cross section) can be written in the following general form 

F2{x, Q^) = — ^ / ^f^^^x - y, z)N{x, y, y) . (2) 
A-K-^aem J 47r 

Here, x and y are the transverse coordinates of the quark and the antiquark in the pair, z is the fraction 
of the pairs longitudinal momentum carried by the quark and y is the rapidity of the slowest particle in the 
pair. Also, $(a; — y, z) is the square of the "wave function" of the photon — the probability that the virtual 
photon fiuctuatcs into the pair with given coordinates and momenta — and J(P{x + y)N{x, y, y) is the cross 
section for the scattering of the pair. This formula is best illustrated diagrammatically through 



\y\ ■ (3) 




The layout is chosen such as to refiect the predominant space time geometry of the collision in the kinematical 
situation at small x, with the a;"*" axis going from lower left to upper right and x~ from the lower right to 
upper left. 

The wave function <f> is well known. It is given for example in Ref. but its explicit form will not be of 
interest here. The focus of the discussion here will be on the "dipole" scattering cross section J N and its 
generalizations. If the quark-antiquark pair is energetic enough the scattering cross section is cikonal and 
determined via the scattering amplitude]^ 

N{x,y)=tr{l-U{x)UHy))^, (4) 



■^Notc that the average over A, denoted by (. . .)a, implies the inclusion of higher order gluonic corrections that are expected 
to be important at small x. It is their influence one tries to trace in the small x evolution equation to be derived below. 



4 



(y dependence suppressed) where U ([/^) is the eikonal phase for the scattering of the energetic quark 
(antiquark) 



+00 

„- A + 



U{x+ = Q,x) = Vcyij>y-ig y da;"A+(a;+ = 0, a;, a;")J (5) 

—00 

with the vector potential in the fundamental representation. One. therefore, has to calculate the average of 
(U[x'^ = 0, x)U^{y^ = 0, y))^ over the hadronic wave function as indicated by (. . . )a- In the present frame, 
the quark and the antiquark move with velocity of light in the negative X3 direction. All the fields in Eq. , 
therefore, have vanishing a;+ coordinate. This will also be the case for all the fields in the rest of this section. 
For simplicity, the a;"*" coordinate will be suppressed in the following. 

In fact this "eikonalization" of operators related to the probe at small x is completely generic, and so, in order 
to describe correlators appearing in any expectation value with the above small x characteristic features, it 
is sufficient to know the corresponding correlations of i/^^^ operators. 

So, to abstract from the concrete example of DIS, one in general is trying to calculate expectation values of 
the form 

(C/(t)®...®C/(tj)^ (6) 

with C/*-^' in the fundamental representation.^ To picture this, just imagine a probe like in Eq. (^) that has 
split into more than just a single qq pair. This again is to be averaged over the field in the wave function 
of the target as indicated by the average over relevant configurations (. . .)a- Clearly one can not expect to 
carry out such averages explicitly without solving QCD, so the best to hope for is to being able to calculate 
perturbative corrections to expectation values like the above. Being interested in small x then means one 
would like to find ln(l/a;) corrections to scattering events. 

This can be done in two complementary ways (or of course any suitable combination thereof). One may either 
try to directly calculate corrections to the target's characteristic field A (for each member of the ensemble 
separately) and this way find the change to Eq. ^ at small x, or one may try to find how the expectation 
values change by trying to modify the wave function of the projectile, which in this setting is the ?7*-^^ content 
thereof. The former approach leads to the JKLW equation, the latter to 1^ and is also the route that I will 
follow here. 

To be a bit more precise, there is a difference in the two approaches concerning the technical setup which, 
although it seems innocent enough at first sight, will lead to differences in the type of nonlinearities one 
is able to include in the end. That difference, although not inherent in writing down Eq. ^ nevertheless 
appears quite naturally and has in fact been used in Refs. |5ff7|,p^. Let me first describe the situation and 
then clarify why it is fairly hard -on purely technical grounds- to step beyond these limitations. 

Any attempt to write down evolution purely in terms of C/^^' correlators appears most natural if combined 
with an asymmetric picture of the evolution: One allows that the target's field, which is responsible for the 
nontrivial exponent of the link factors may be large. At the same time the field generated by the projectile 
is assumed to be small, it plays no role in finding the U^^"^ correlators. Clearly such an assumption would 
break down if evolution were to generate large fields on the projectile side. As it will turn out gluons 
are produced well spread out over rapidity and hence the projectile's field may be considered to be small 
throughout. There is another problem, however, that will become crucial as projectile rapidities approach 
target rapidities and the two objects -components in the tails of their wavefunctions really- cease to be 
distinguishable. Then this picture will clearly break down. Evolution in the JKLW sense, on the other hand, 
evolves the large field of the target into an even larger field, taking all the nonlinearities into account already 
in the emission process. Technically speaking the fluctuations propagator around the dominant background 
flcld configurations contains additional nonlinearities that go beyond the "independent emission" scenario of 
the weak field perspective imposed on the projectile side. 

^This is completely general as one may write any adjoint links as a combination of two fundamental ones by virtue of 

Uab = 2tv[rut>'u^. 
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To overcome this limitation in the projectile evolution picture, one would have to allow both fields -the target's 
and the projectile's- to be large. Due to the nonlinear nature of QCD this would require to determine the 
relevant dominant field configurations due to two strong sources and find the fiuctuations propagator around 
such configurations, a feat of which not even the first step has been accomplished analytically. 

From this perspective, there are three degrees of complications with correspondingly increased amounts of 
nonlinearities included: Evolution from a projectile with a weak field in the presence of a strong target field 
as done in [p|-|7|,[l^; evolution from the target with a strong field as in the JKLW case and evolution with 
both fields strong, which has not yet been looked at in earnest. 

Returning from this general discussion of the scope and content of different approaches to matters at hand, 
let me continue to set up the calculation in the projectile evolution perspective. Here the idea is to calculate 
small X corrections to correlators of the type Eq. (^) and infer their change from the resulting equations, 
without explicitly tracing how the target fields change as one proceeds. 

As one is interested in the behavior of arbitrary correlators of the form Eq. (^) it would be most economical 
to find the generic behavior of the generating functional of such objects 

Here 

Sll,[A.J\J] = jd^x[tr{{JlfU4A+])+tr{JiUi[A+])] (8) 

is an external source term,^ correlators (Eq. (||)) are extracted via J^^^ derivatives and the {. ■ ■)a average is 
taken with the QCD action. 

So far the physics of small x is encoded in the type of correlation functions considered -exclusively correlators 
of link operators ?7^^'- all the rest is mathematical convenience. The sole purpose of this machinery is to 
get a clearer, more streamlined picture of what one can do to calculate these correlators and how to do it in 
practice. Simple, diagrammatic interpretations are close at hand. 

As it is obviously impossible to perform the field average in one go and the form of iSg^t [A, , J] has the 
small X limit already built in, one will look at Eq. expanding the gluon field A around h (whose correlators 
are assumed to be known) with fluctuations only taken into account to order as ln(l/a;). That is to say that 
one will expand around h to one loop accuracy and select the terms carrying a ln(l/2;) factor. This way one 
will be able to infer the change of correlation functions as one lowers x. This is the second statement about 
physics or rather what one can learn about it through such an approach. 

Now turn back to the actual calculation. At one loop one needs at most second order in fluctuations:^ 

1 + 5A^^ + hA^±SAy4- + 0{5A')) e^'!^^^^' '-'^sa 
obx 2 obx oby J 

+ LsA^1^6Ay^e^'i^^''''^%.sA + ... 
2 obx oby 



l^e&^J'-\ + ]^{{5A^5Ay)sA[h] 



(9) 



X 



2-^5^|j6, Jt, J]^Sll,% J\ J] + J] 

X y X y 



Note that I purposely write tvA^B to unclutter notation in what follows, not trAB, having in mind that ,4 . . tvA^B 



S 

6A., 



- as discussed above 
In writing Eq. (^) one 
propagator in this background field and has been used repeatedly 



^ In writing Eq. (|^) one has anticipated that (^aZ) g^[b] = as iji^ |h^ free case. This is a consequence of the structure of the 
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The second term is the calculable perturbative correction to Z[j\ J] the generating functional for all gener- 
alized distribution functions of the target, an object which can not be calculated as such with present tools. 
Arbitrary correlators of course are extracted according to 



5s,iAu,u\j\j] ^ f^te5f,[c/.c/t,jt,j] ^ ^^^^^ 
J 

5J'^ 



or 



_5_ 5 5 5 



A^TTt • • • TTF^t^^' J] = {Ul ® ■ ■ - ^Ul^^U^^^ . . .®U^)^ . (10c) 



To understand the individual terms in Eq. (||) in detail requires the calculation of the one loop corrections to 
path ordered exponentials C/[6]^t) the presence of a background field b of the form Eq. (|l]). Before diving 
into this it may be helpful to specialize once more to the DIS example from above to illustrate the physics 
content of the terms in Eq. (Q). Here one needs to look atR 



= (tr(l-[/46][/t[fo]))^ 



5 5 

quantum corrections 



The averaging procedure (. . . )t, is a; dependent as the average over SA brings In l/x dependent terms denoted 
by "quantum corrections", a special case of the second term in Eq. (|^). This sets the task: // one is able 
to calculate these quantum corrections before taking the b average, that is, for all relevant b, one can deduce 
how the [/(t) correlation functions or, as will become clear in Sec. the weight Z[U, C/t] ffiaf; defines {■ ■ ■)b 
evolves with x. 

Clearly there are two generic types of corrections corresponding to the two terms in 

\{5AJAy)sA[b] ((^-^flAb, J\ J\ J] + ^^^/o^[b, J\ J 

X e^-[^-^'-^l . (12) 

The first term represents contributions where the gluon propagator connects two quarks [U s) or antiquarks 
(C/ts) as well as a quark to an antiquark. In addition there arc pure self energy corrections dressing one quark 
or antiquark line instead of connecting two of them, represented by the second term. 

Particular correlators are again selected taking any number of J derivatives and setting J to zero. One point 
functions get contributions from the second term only, while both terms contribute to anything with more 
than two J^^) derivatives. Clearly it is sufficient to calculate 



1 5 



2 wit) W^t) 



j(t)=o 



\{5A^5A,)sA[b] J\ J]j^/tlt[b, J\J]] (13) 



and 



5J., 



(t) 



,/<t)EEO 



\{5A^5A.,)sA[b] (-^^/Zh J\ J] ) (14) 



to completely reconstruct Eq. ( |l2| ) as these terms are precisely second respectively first order in j't). Eq. (^2|) 
then defines the change in all other correlators and one has reached the goal of finding the fluctuation induced 
corrections. 



^ After integration over impact parameter (simply x + y, if one puts the target at the origin in transverse space) this is the 
dipole cross section that features prominently in calculations of 7*p or 7*^ cross sections at small x [|7|,|l9|,|20|. 
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The task is clear and the terms appearing in the actual calculation are best visualized diagrammatically. For 
the gluon exchange diagrams corresponding to Eq. (Ol) one defines 




(15) 



As already seen in these split naturally into x~ ordered contributions when one combines the structure 

of the vertices (the b derivatives of ^cxt in Eqs. (|l^) and (^4|)) and the gluon propagator in the background 
field. Take x''^ a-s an example:^ 




"a ln(l/a::) cr^ 




(17b) 



Again, the various components of x correspond to gluon exchange diagrams while the as are self energy 
corrections. As for the time ordered ones, there are those with and without interaction with the target 
(with the gluon line crossing the axis or not respectively as shown in Eqs. (|l^ ) and (17)). These latter 
versions determine the J/^^^ dependence of x and tr in detail. They have been calculated several times using 
different methods. The common feature of all of them is that diagrams without target interaction leave the 
number of [/'^'s invariant while those with target interaction insert an additional adjoint C/'t) at the point of 
interaction. This implies a nonlinearity in the evolution: Any correlator with a fixed number of fundamental 
representation U^^^s will couple to other correlators with two additional factors of fundamental representation 
[/'t)s in each step in x. 

It is probably most efficient to read off the results for x and a from Ref. II or ll8|. For convenience, define 
the integral kernel 



{x - z).{z - y) 
(x - zf(y - 2)2 



(18) 



'^Representations for Xxy Xxy ^.nd Xxy result from reversing the quark lines aecordingly. 
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This allows to write explicit expressions for the components of x 



o^sbcVyhki := ^ jd z IC^^y (19a) 
X {[U,UlUUlU^]kj + [U^Ul]u[UlU,]kj - [U^Ul]u6kj - Su[UlU^]kj 
[xllh ki := as [Xyx] w 'J = ^ / ^y^"" (1^^) 



asx'''^{x,y)tjki ■■= -77^ jd^zJCx^y (19c) 



OisixVyhki := jd^z JC^^y (19d) 

as well as those of ct:^ 

- - ^ /'^'^ /C,.x([C/,]ytr(C/,C/t) - iV,[C/,],j) (20a) 

--^-^ Jd'^ ic^.4PlhMulu.) - NcPlh) • (20b) 

At this point all the calculational effort pays off: Using the above definitions for x ^md ct, and remembering 
that S^^^[A, j\ J] depends on A via [Z*-^^ only, recasts Eq. (^) as 

+ \{{5A,5Ay)sA[b] {^^/lUb, J\ J]^flt[b. J\ J] + -^^^/e!t[b, J\ J] 



a. 



(21) 



The logarithm comes from an integration over 5 A with momenta in a finite interval of x values. Accordingly 
the equation has the form of a finite difference equation for Z[j\ J] with respect to Inl/x. 

This gives rise to the final RG equation in differential form, which I write using condensed notation with 
J=[J\J) and^ = (^,^)as 

-Z[J] = as \ -JuJvXuv[T-f] + Ju^u[t^] \ Z[J] (22) 



Inl/x [_2 oJ oJ 

where u and v now also stand for q and q in addition to color and transverse coordinates. Eq. ( |2^ ) represents 
a full set of evolution equations for all correlators of C/'^'^^ that can be extracted by j'^) derivatives evaluated 
at J^t-' = as demonstrated in Eq. (|ll|) for the dipole cross section.^ It is probably worthwhile to recapitulate 
the surprisingly few steps involved in arriving at the quite complex result displayed in Eq. (p2[). 

• The first step consists of writing down the generating functional of all the correlators important at 
small X. It is then a small step to recognize that quantum corrections at small x can be calculated 
directly as corrections to this object rather than for all the correlators separately. 



®Note the minus signs in K^zx = — l/Cf — ■ 

^The convention for index contraction again is defined in order to simplify notation. Products are defined as in jto-' 
[jt] [g-'Jij with integration over transverse coordinates impfied 
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• Simply expanding in fluctuations allows to identify the main players in the evolution x and a. Their 
calculation requires some effort, but also that can be streamlined considerably | [T^ . 

• Even before calculating their values, with just their basic presence known, one then may write the RG 
equation in terms of x a-nd a. 



The physical content of the RG is equally transparent. Eq. ( |22[ ) describes an infinite set of coupled RG 
equations for correlators of U fields averaged over the dominant color fields in the target G^^[b]{w) 
b~^{w) = P{w)S{w~). As is clear from Eq. ( pi|) one can describe this in two complementary ways: First as an 
evolution in which the target is left unchanged -the (. . .)b average stays the same in an evolution step- while 
the corrections induce correlations between different operators characterizing the projectile wave function.]^ 
This description has been used many times in the literature. Second, and this interpretation is closer to the 
spirit of Eq. (|2^), as an evolution in which the correlators of U operators are changed through an infinite 
set of coupled RG equations. As one is talking about correlators of U averaged according to (. . .)b this 
implies a change in this averaging procedure as x runs. In turn this is nothing but a change of the dominant 
configurations contributing to this average and hence a change in the target wave function. 

Any way one looks at it, the only operators involved are correlators ofU — U[b] averaged over b. Accordingly, 
during the entire evolution, one will encounter nontrivial correlators only within the targets -x dependent- 
radius. 

During the course of this work it will become possible to make much farther reaching statements about this 
evolution. For now I will return from generalities to concrete examples that may help compare individual 
expressions to Ref. and |^] and should give a better idea of the structures involved when trying to represent 
Eq. (E3) as a system of coupled equations. 



The one point functions evolve according to 



d 



dlnl/x SJ, 



d 



91nl/a; (5 J; 
or equivalently, in the language of correlators 

d 

d\n\/x 
d 



as all—] Z[J] 
as aii^] Z[J] 



as{a^JU])i, 
as{ai[U])t . 



(23a) 
(23b) 

(23c) 
(23d) 



dlnl/x 

To emphasize this point once more: since the as are nonlinear in C/^^^ these equations do not close among 
themselves but couple (by induction) to all other evolution equations for higher n-point correlators. 

• An example for the evolution of a two point correlator (there are 3 more generic cases involving the 
other components of x) is given by 



d 



d\nl/x SJl SJy 



(24a) 



j=o 



5 J SJii 



or, again equivalently, 



d 



:(C/x ® Ul)t = as {x%[U] + al[U] ®Ul + Ux® al[U])^ 



91n 1/x ^""^ "^'^ \^xyi^ J ' "xi-^ 1 ^ ■ -^--"T/L-j/f,- (24b) 

This last equation contains the evolution of the dipole cross section alluded to in Eq. (11), a connection 
that will be explored in somewhat closer detail in Sec. ^ 



^"This can not be true down to arbitrarily small x where projectile and target strongly overlap and limits the region of 
applicability of this particular RG equation, see jl^ . 
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3 Limiting cases 



The full RG equation is quite general and contains various limiting cases. Two distinct types are of particular 
interest, 

• the first one which still captures saturation and unitarization effects, the dipole evolution for 'y*A 
collisions based on a large Nc argument and 

• a second regime in which expansions at small "densities" simplify the equations. This may be called 
the BFKL limit. For obvious reasons this regime will only be useful to gain insight at moderately small 
X before saturation effects set in. 



To turn to the nonlinear A^^ dipole evolution for "f*A collisions first, recall Eq. ( p4q ). ft allows to extract the 
equation for the dipole cross section without much effort. One reads off that 

X 2 [tr{UlUMUzUl) - NMUlUx)] 

while the sigma contributions are given by 

H]., H],. =^ /^'" (^^f)^(Sr^ HuluMUyUl) NMulUy)] (25b) 

H],. [Uy],, J'^' ^z-lnx-zY HulUyMUlU.) NMUlUy)] . (25c) 

Adding these according to Eq. ( ^4b| ) leads top| 



d 



— {tviU^UD), = a,tr{x%[U]+al[U]Ul + U^al[U])^ 



:{tr{UlUMU.Ul)-NMUlUx)),. (25d) 



dlnl/x 

27r2 7"" [x- zYiz-yf' 

This equation exhibits the well known dipole kernel and turns into known closed evolution equations, if 
one 

• factorizes at large according to {tr(UlUx)tT:{UzU^)) ^ N^^oa^ (i^ciUlUx)) ^(^tr{UzU^)) ^ to decouple 
the equation for the two point correlator from higher order ones to close the equation. 

• expresses the resulting equation in terms of the dipole scattering amplitude N^y '■= (tr(l— ?7xC/^))^/Afc) . 

As a result the RG equation of Ref . ^ for DIS emerges 

3 N f (x — 1/)'^ r 'I 

m^x'^-y - -if J'' \J.)4-yA ^''- + ' ^^^^ - ■ ^''^ 

This equation has hallmarks of an equation that incorporates saturation effects. Schematically this is quite 
simple: Starting from small N, the linear terms are the only relevant ones and grow according to the BFKL 
kernel they are subject to. At the same time the N'^ term grows as well and will catch up with the linear one, 
at which point one has reached saturation. Although transverse dynamics will be important for the details 
of when this happens, it will not alter this basic mechanism. 

Note that theLcancellations leading to Eq. (25d) provide a cross check on the relative sign of the x'^'> and both the <t 
contributions via ItI. Unfortunately the relative signs of the other x components are not fixed by this, but there will be a 
surprise waiting when I come to look at the force term in the FP formulation. 
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The second limit of interest is the small "density" limit in which one may start expanding all U^^^ factors 
around 1. Alternatively one may attempt to decouple the RG equation hierarchy by expanding in powers of 
U^^\ Both of these perspectives will have only a very limited region of applicability as will be seen below 
and will most likely yield similar results there. 1 will limit myself to the latter. 

To see how these expansions affect the RG, define the generating functional for connected correlators via 

Z[J]=e-^[-^l (27) 

and view G[J] as a power series in J = ( J^, J). Similarly think of x[?7, C/^] and ct[{7] as a series inU — {U, W). 
For example for x one writes the expansion 

Xuv[U] = J2 ^xl"L,...a„t/o, ® . . . ® f/a„ (28) 

with Xuv:ai...a„ thc n-th derivative of Xud[C^] with respect to its argument and similarly for a and Q. 

Truncating x and a at lowest order (quadratic and linear in U respectively) they reduce to real and virtual 
parts of the BFKL kernel in the "dipole" formulation. One then finds 



d 



d\nl/x 



g[j] = as |^J«J.xi'J;„;3( - g-MJ] + 5-AJ] g-AJ]) - Jval^lGAJ]] 



(29) 



This still represents a whole hierarchy of equations; one equation for each power in J. Assuming a vanishing 
one point function, the hierarchy starts with the BFKL equation for 5^^^. The equation for the three and 
four point functions then get feedback from the BFKL equation via fj'^', a pattern of lower order n-point 
functions coupling into the equations of higher order ones that continues on to arbitrary Q^"\ This is the 
very same pattern that shows up in [|2|, although here one is looking at "projectile" evolution. Had one 
started this type of truncation scheme with the JKLW equation, one would have directly reproduced the 
"target" setup. 

This identification together with what has been established in Rcf. should now cover all presently known 
equations for small x evolution to leading order in a^lnl/x in which target and projectile are treated 
asymmetrically. 

Expanding the t/'^-' factors around 1 will lead to very similar equations for the lower order correlators, 
although clearly higher orders will show differences for the two approaches. However, it is probably pointless 
to try and contrast the two versions based on these higher orders. Clearly as it becomes important to include 
such contributions there will no longer be any justification to truncate the full evolution Eq. ( p2|) in any form 
and both versions lose their validity, although one might outlast the other by some margin. This -as well 
as potential saturation properties of these truncations- is not of much interest in the present context, as the 
new results to be developed will be for Eq. ( p2[ ) directly, without the need for any truncation. For this same 
reason I will not comment on other cases, such as the doubly logarithmic limit, but rather refer to Rcf. 
for the JKLW case and Rcf. for a comparison between the target and projectile evolution pictures also in 
that limit. The result there nicely mirrors the the hierarchy of increasing complexity described at the outset 
of this derivation of evolution equations. 



4 A Fokker- Planck form for the evolution equations: Probability 
distributions &; fixed point equations 

If one recalls the definition of Z[J] := (e'^''t[''''^l)(, and the initial idea of interpreting (. . .)6 as a statistical 
average with associated statistical weight Z according to 

{...)t^{...)u = jD[U]{Z[U]...), (30) 
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one may try to rewrite Eq. ( |2^ ) as an equation for Z. Indeed -remembering that 5^xt[^: "^1 — '^cxtl^i "^1 
actually given as a functional of U'^'^^- 



d 



dhil/x 



D[U]Z[U]e 



, I D[U]Z[U] 



(31) 



Using the fact that if one treats U and W as independent variables one finds, conjugate to Eqs. ([l0|), that 

(32a) 
(32b) 



_6_ 

SU 
6 



Hence, Eq. (|3l| ) can be rewritten as an equation for Z. Note that now both t/'^t) ^ GL{Nc,C), not SU{Nc), 
otherwise the condition UU'' = 1 would invalidate Eqs. (^2|). The restriction to SU{Nc), which of course has 
to be present for consistency of the whole approach, will have to come from the properties of x f-nd a and 
hence the nature of evolution for Z. Leaving this for later, the evolution equation for Z is given by 



d 



91nl/a; 



Z[U] = as 



1 S 



2 SUu SU^ 



Xuv[U]-—au[U]\z[U] , 



(33) 



where the J/*-^^ derivatives act on both the kernels x,a and Z. This equation exhibits exactly the same 
structures as the JKLW equation, although the variables and details of the kernels are, of course, different. 
Note that this has the form of a Fokker-Planck (FP) equation, simply due to the fact that the original 
equation as represented in Eq. (31) had only terms quadratic and linear in J. 



As with any other FP equation, the fact that the right hand side is a total derivative implies that the 
normalization of Z -its U, integral- is conserved under x evolution and, if one chooses this norm to be 1 , 
allows to interpret Z as a probability distribution at any stage of the evolution. 

Clearly this formulation is completely equivalent to the original version involving the generating functional 
Z. Instead of obtaining equations for correlators by functional differentiation, one now extracts them, again 
as in any FP formulation, by multiplying with the desired mononomial in t/'^-' followed by integrating the 
result over 

Let me conclude this section by rewriting the FP equation Eq. (^ ) in what I would like to call its canonical 
form by pulling one of the derivatives in the x term to the right and factor the other derivative out to the 
left. The FP equation then reads 



9 In 1 /x 



Z[U] 



1- rrnJ- 



SU,' 



Xuv[U] ]~au[U] 



Z[U] . 



(34) 



Technically, this form is particularly useful when searching for fixed points, but the reason for this is of 
course rooted in physics: Written in this form the terms have the explicit interpretation of a stochastic 
and an external force term. Only the latter would be capable of balancing the diffusion caused by the first 
term and could lead to fixed point structure that does not correspond to pure diffusion. However, upon 
translation of Eq. (|3^) into a Langevin equation, it is still a that that shows up as a (possibly "apparent" ) 
force, while x determines the correlator of the purely stochastic fluctuations. However, I will delay writing 
down this Langevin equation -the prerequisite for any numerical simulation for such evolution- until a better 
understanding of the properties of the RG equation has been reached. 

Any reader familiar with typical examples of FP equations, be it in classical physics, stochastic quantization 
and its applications to numerical simulation of field theories or any others, will by now have begun to wonder 
about the the "equilibration" behavior of the FP equation (^3|). In the present RG context, this is nothing 
but the question for the existence and nature of the small x fixed point of the evolution. 
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A sufficient condition for the existence of a fixed point of the RG is simply given by the vanishing of the 
driving term of the RG, i.e. the r.h.s. of Eq. (pj). In other words, a fixed point exists if 



Z[U] (35) 



has a solution. Implicitly this determines the fixed point value of In Z via a functional differential equation in 
which one has to integrate back from the force term to get the "potential" that determines the "equilibrium" 
or fixed point. It depends on the nature of this solution whether the fixed point in question is attractive and 
hence reached by evolving any arbitrary initial condition or not. Besides the mere existence of a fixed point 
it is this fact which determines whether Eq. (^4|) predicts universality in the sense of target independence of 
QCD in scattering events at small x or not. Instead of pronounced regularities, a repulsive fixed point would 
imply runaway solutions and indicate an increasing target dependence of measurable quantities with smaller 
and smaller x. That way the force term on the r.h.s. of these equations will decide the general nature of 
QCD cross section at small x in the leading Inl/x approximation employed. 

To answer any of these open question one will therefore have to carefully study the properties of both the 
force and stochastic terms. 



5 A case of Brownian motion? 

Indeed, there remains quite a bit to be learned about the character of the evolution described by the above 
FP equation before one can satisfactorily answer this question, but a first observation is very easily extracted, 
just by looking at the force term. 

To this end one has to evaluate U'-^^ derivatives acting on x and, quite surprisingly, one finds 

U4rXuv[U]) -au[U]=0 . (36) 



2\SUv 

(For both the q and q components.) During this calculation, which is exemplified in App. one observes very 
intricate cancellations of infinities between the derivative terms before the finite remainder cancels against a. 

As a result, the force terms vanish completely and one is dealing with a Fokker-Planck equation with a 
nontrivial kernel, but without a force term. The cancellation of the latter is not at all an accident. In fact 
it relies on a very deep relation between the x ^'Hd a terms. Tracing back what is left of this cancellation 
in the BFKL limit, one finds that there it is responsible for the infrared finitencss of the BFKL kernel. 
There gauge invariance guarantees that one finds a cancellation of "real" against "virtual" diagrams in the 
infrared. The situation and the cancellation just observed here is even more general, as it goes beyond the 
BFKL limit, that is the limit of the RG equations where U is in the vicinity of 1 to include the full set of 
nonlinearitics of this evolution equation. 

With the force term vanishing, the Fokker-Planck equation now deceptively looks like a case of Brownian 
motion: 

^ -Z[U]=a,lJ-XuAU]JrZ[U] ■ (37) 



dlnl/x ' ' ''2 5Uu SU, 

One might be tempted to jump to the conclusion that Z ~ const is the "equilibrium" limit in standard 
Fokker-Planck parlance - the fixed point at small x in this case. In particular as one is forced to take U 
and derivatives independently, thus violating for instance the condition that UU^ = 1, one seems to be 
in danger of losing all memory of the fact that QCD has a SU (Nc) gauge symmetry, a possibility strikingly 
at odds with what has just been identified as the source of the cancellations leading into this very dilemma. 

Barring any hope for an inconsistency an unfriendly reader might harbor, there must be further structure 
hidden in Eq. (|3^) that prevents this from happening and the only place left for that to hide is within the 
nontrivial kernel x- 

^■^This statement is strictly true only in the "target" formulation. For the present "projectile" or "dipole" formulation the 
same statement refers to a cancellation of unphysical ultraviolet divergences. 
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6 Uncovering the gauge group 



6.1 Independent derivatives? RG flow and Lie derivatives 

The fact that in writing the Fokker-Planck equation one had to take independent derivatives with respect 
to U and W , raises the question of how the evolution knows about the phase factors being SU{Nc) and not 
some arbitrary complex Nc x Nc matrix. Generically there are two options: 

• To outline the first possibility let me turn to simple mechanics analogue, a Fokker-Planck/Langevin 
problem in a 2d harmonic oscillator with a stochastic term acting only in a Id subspace. This in 
turn implies that the system is completely deterministic in the other direction and the restoring force 
there will asymptotically move the particle into the line with stochastic perturbations, the rate of 
approach being proportional to the force. The probability distribution (for all times) will factorize into 
a deterministic and a truly probabilistic part. The deterministic part has the form of a time dependent 
delta function following the particles trajectory from its initial value into the minimum. If the force 
becomes infinite (say by making the potential infinitely steep in this direction) this factor will become 
time independent and fix the position at the minimum in the non-stochastic direction. This is to say 
that an infinite force signals a constraint in the system. 

Although this looked like an attractive option when discovering the infinite contributions to the indi- 
vidual terms in Eq. ( |3^ ) it is emphatically ruled out as a consequence of the remarkable cancellations 
occurring when adding all of them together. 

• The second option is simpler and more elegant: Imagine the forces (both deterministic and stochastic 
alike) to act only "along" SU{Nc), not out of it in the sense that if you start with an initial condition 
"within," evolution will not push the system "outside." Clearly this would have to result from intrinsic 
properties of both x ^-nd a, although in the present case the burden lies solely with x &s the only 
independent ingredient left in the equation. 

In this scenario evolution would separate the whole (apparent) configuration space into invariant orbits. 
Given an initial condition in one of them, evolution would keep the system within this same orbit. To 
achieve what is needed for consistency, one of these orbits would have to be SU{Nc) as embedded into 
the complex Nc x Nc matrices or GL{Nc, C). 

Given this property, one would be able to reduce the above equations to the orbit of the Fokker-Planck 
operator on an invariant subspace. In other words it should be possible to write the FP equation in a 
more compact form within what is of interest here: the physical subspace of SU{Nc) matrices. This 
will be done below. 



The first clue that a separation into invariant orbits, with SU (Nc) being one of them, might indeed be realized 
is the observation that 

Xuv[U]-^ [U^UlUn ^ =0 (38) 
oUv c/c/t=i 

for both q and q components. 

What one needs in fact is, that for physical distribution^^ of the form 

Zphys[U, C/^] = (5(C/C/t - l),5(det U - l)Z[U] , (39) 

the constraint factor S{UU'^ — l)S{detU — 1) is left invariant by the evolution operator, i.e. that in some 
sense one may commute it though the delta functions. 

As it turns out the evolution operator can be pulled through the constraint and does simplify tremendously 
during this procedure. Only then will one be able to see if the Brownian motion concept is actually realized. 
Pulling the evolution operator through the constraint is in fact the main technical problem left. 

^''all others would be illegal and spurious in any case 
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In order to get a feeling for what is involved it is useful to collect some ideas of what to expect if this were 



actually possible. The main step in rewriting the evolution within the group consists of turning -Aj and ^ 



from independent to dependent derivatives via UU^ = 1. What is to replace them? 

Pondering this question a little (c.f. App. |^ for a pedestrian approach), one inevitably ends up considering 
the natural derivatives along SU{Nc), its Lie derivatives. They would be given by^ 

^u^=Pn4^^-[-t^U-%j^^^, (40) 

so that as a first step in this direction one would like to reexpress things via [Uf^jij-^-j and [—t^U^ijj^.. 
instead of -^^j and -^jj^y Indeed, if one now just tries out what this transformation will do, one finds 
remarkable results for the transformed components of x- 

To achieve this technically, one needs to be able to represent all of 5/5u\J^s components. Accordingly one 

has to replace the index a e {1, . . . , iV^ - 1} by A e {0, . . . , TV^ ~ 1} and supply f = lJ_s/2Nc 1, such that 
tr(t'^<^) = 5'^^ /2 for all A, B. Then the transformation is easily accomplished (see App. |c|for details). 

As already advertised, the results arc quite striking: 

• in SU{Nc) all four of the transformed components are equal. 

• in SU(Nc) all of them vanish if cither A or B is 0. Only octet components survive. All matrices involved 
now are explicitly in the adjoint representation. 



• The unique common form even factorizes naturally into the "square" of a much simpler factor. The 

[ 



result will be denoted by X™[?7] and reads 



The diagonalizing factors of x deserve a name for further reference 

(42) 

If the FP equation can be written entirely within SU (A"c) ~and this will be demonstrated below- one may 
also expect that a description in terms of a Langevin process with a completely decorrelated Gaussian random 
noise becomes possible due to the factorized form of (^ij) via (^^. 

6.2 A Fokker Planck equation within the gauge group: The small x fixed point 

The results of the above transformation of x are so suggestive that it seems obvious what the form of the FP 
operator on physical configurations -modulo some clairvoyance as regards the prefactor- has to be: 

as^iVtj^rJ'y^'^lj, ■ (43) 

The chain of argument to prove this is the following. One sets out to show that 

^ Xuv4r ^iuu^ - i)'5(dct[/- 1) = 5([/c/t - i)5(dct;7- 1) ^irl'^l,, ■ (44) 



This amounts to verifying that 



-"^^It is easy to see that this is the correct form as both variants amongst themselves and with each other satisfy the SU{Nc) 
commutation relations at each point. 
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for arbitrary mononomials U^l (8) . . . (E" U^i , a task that once more is best performed via a generating 
functional. This strategy has the added benefit that it makes it readily apparent how the two operators 
lead to the same hierarchy of equations for the correlators which provided the starting point for all these 
deliberations. The details of this exercise arc given in App. |^. 

This leads to the central result of this paper, the RG/FP equation on the physical configuration space, where 

Zphysp, C/t] = (5(C/C/^ - l)(5(dct U - l)Z[U] . (46) 

The constraint factor may now be absorbed into the measure to form a functional Haar measure according 
to 



r(t) 



D[U] := D\U]5{UU'^ - l)5(dct U - \) 



(47) 



This constitutes the natural measure to use for averaging correlators of the form (Uw^ <?) ■ ■ .® Uw„ )jj with 
the probability distribution Z[U]. The RG/FP equation simplifies to 



Z[C/](x) = a. l^y'lIJly^<, mix) 



(48) 



This formulation is fully equivalent to the initial one, but with a large amount of redundancy removed. 

Contrary to Eq. (|3^) the evolution now docs have the interpretation of Brownian motion in the ordinary 
sense: 

To see this write the evolution equation in terms of what is canonically called the Fokker-Planck Hamiltonian 
Hpp as 



where 



Now dcfine|^ 



to recast Hpp as 



d 



dlnl/x 



Z 



-HfpZ 



H 



FP 



pa 



ab 



idi 



2-K 

92 



xy 



(49) 



(50) 



(51) 



H 



FP 



47r2 



ptp 



(52) 



Hence, the spectrum of of iJpp is bounded from bclow.^ This is what one was hoping for: Now the system 
will "equilibrate" , the small x fixed point will be attractive: A finite evolution step is described by 



Z{V) = e-^-(«-^°)Z(yo) , 



(53) 



leading to a damping of all but the zero eigenmode, which will become the fixed point. This fixed point 
however is trivial: With iJpp acting as a derivative, the asymptotic solution to this evolution equation is 
simply 



Z[U] 1 , 



(54) 



^^Note how crucial the sign here becomes: If it were negative, the Langcvin formulation would break entirely. Instead of getting 
random kicks in the Lie algebra, one would see an additional factor of i there and such a formulation becomes meaningless. It 
is instructive to note that the wrong sign here would simultaneously imply a decreasing BFKL pomeron as one easily discovers 
when tracing this sign back to what it corresponds to in Eq. (|26[). A growing BFKL pomeron and positivity of -ffpp are the 
same thing. 
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where it has been assumed that the Haar measure has been normahzed to yield unit probabiHty within 
SU{Nc)- While Eq. (jsj) obviously is a zero mode of Hpp, it is much less obvious that it is a unique one. 
The argument here rests on the fact that the Lie derivatives in Hpp induce local gauge transformations with 
respect to the transverse coordinates, while any change in the ± plane is immediately related to the coarse 
graining step of the RG procedure. For the details, see App. ^ 

It is definitely worth taking a step back to contemplate this result; The fixed point probability distribution has 
no U dependence. Its physics content lies in the fact that the distribution of ?7's completely fills configuration 
space and the "phases" randomize entirely, a standard Brownian motion scenario. This may sound trivial, 
but in fact it will have intriguing consequences. On the one hand this is the root of any universality argument 
i.e. the hope for target independent features at small x. On the other hand, as will be shown in Sec. |^, this 
limit of infinite energy is not the most interesting piece of information to be learned about the asymptotics 
of this evolution. Already the fact that this limiting solution exists and is attractive leads to regular features 
in a larger region of energy, well before the fixed point itself is reached. To expose this additional structure, 
however, one needs to establish a few further insights. 

One is, however, already able to get a glimpse on how the damping exhibited by Eq. ( |53|) can be reconciled 
with the known strong growth of the BFKL limit. This in fact will expose the main unitarization mechanism. 
Note that it is the dipole kernel K which governs t he B FKL limit of Eqs. ( ^5d| ), ( p6|) not JC, the integral 



kernel in x. They are related according to (c.f. Eq. (25a)) 



JCxzy — Z 



1 r ~ 



2 



(55) 



In contrast to /C, K is not a square and in fact not bounded at zero. This is in direct correspondence to the 
same lack of boundedness at zero of the BFKL kernel. This observation sheds some light on the unitarization 
mechanism and the lack thereof in a strict BFKL truncation. To see how this comes about take Eq. ( p6| ) , but 
before factorization. Define the dipole operator A/'xy := tr(l — UxU^)/Nc (not the amplitude). Its evolution 
equation is given by 

and clearly couples to all higher order ones through the nonlinear term. Now think of the full evolution 
equations as an infinite set again, more precisely, think of it as a matrix representation on the space of 
correlators. Eq. (|5^) only reduces to the BFKL equation if one ignores the nonlinearity and hence the way 
this particular piece of the matrix equation couples to the other components of this matrix representation. 
That is what breaks positivity and leads to the runaway solutions. 

This, however, also explains -along the same lines as in the large Nc truncation of how a unitarization 

mechanism might work here: As long as the correlators are close to the free ones, the nonlincaritics are 
unimportant. The growing mode of K can work its magic and violently push the system towards a situation 
in which the nonlinearities become important. Then "positivity" of Hpp takes over and slows down all further 
evolution drastically. 

Complementary to Eq. there is one thing one knows about the initial conditions, and that is that all U 
factors outside the transverse tail of the wave function are unity. There 

Z[U]^6{U-l). (57) 

This and the notion of what is outside the target will change during the course of the evolution, as x has non 



vanishing contributions everywhere, due to its nonlocal nature, c.f. Eqs. (41), (|42|). Together with the initial 
BFKL like growth, this will be an integral part of the way unitarization is reached. I will come back to this 
after having spelled out the Langevin equations. 

Eq. (|^ together with Eq. (|5^) implies a full knowledge of all correlators as one enters the saturation region 
with X —f 0, provided some information about the transition region and hence the transverse tail of the wave 
function is put in via the initial condition. 
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I will come back to use these correlators below when discussing the consequences for unitarization in DIS. 
As seen above, impact parameter space and relative coordinates play a very different role in the evolution. 
The latter only enter in a scale invariant manner into the kernel. The interplay of these two facts will lead 
to scaling laws with respect to target "size." As this applies to the evolution of all correlators, it will also 
have consequences for the DIS cross sections. 



Langevin equations: the basis for numerical simulations, scaling 
laws 



Before discussing the generic form of correlators in the saturation region, there is one more useful perspective 
to be explored that can be harvested without much effort, directly from the the wealth of methods in the 
stochastic physics tool box. Besides providing a generic method for numerical simulation, it will also provide 
a very intuitive physics picture. This will directly lead to scaling laws governing the whole of the evolution 
and provide a very nice, novel perspective to view the unitarization mechanism. 

There is an extensive literature on numerical simulations of field theories with stochastic methods, that can 
be adapted to simulate the generic small x evolution equation Eq. (|4^). This is an enormous step beyond 
what was feasible up to now . The key tool to do so lies in the possibility to describe the same stochastic 
process encoded in Eq. ( ^8| ) via a Langevin equation with Gaussian noise. 

To prepare translation into Langevin form let me write the RG again in standard form 



d 



d\n{l/x) 



Z[U] = -a^Vl \WyXty - {\^'yXty)\ Z[U] 



and define 



2 y-^^v ^ 



^2 7r2 J {x^zY 

(the "~" indicating adjoint matrices and traces). The Langevin equation then reads 
d 



9 In 1 /a; 

where both uj and ^ are Gaussian random variables with correlators determined according to 



(...). = detx'/2|i5H(...)e-5-*-'- and {. . ^ J D[^] {. . .) , 



(58) 



(59) 



(60) 



(61) 



Also note the factor of i which is essential to render this an equation for an infinitesimal change of an element 
of SU{Nc). In fact the components of uj^ can be directly interpreted as the "angles" parametrizing a local 
gauge transformation in transverse space. 

It is of particular importance to note that the possibility to formulate the stochastic term via a completely 

(2) 

decorrelated Gaussian noise that is to say with (C^^^) = S^^Sxy, reduces the numerical cost for a simulation 
like this considerably 

Note how the Langevin equation induces transverse spreading: Even at points initially outside the target, 
where Ux starts out as unity, the noise term generates nontrivial phases due to the nonlocality in [f^^];. This 
growth is of the order of ^ylnl/x according to the typical scaling behavior of all Brownian processes, by a 
textbook argument. How does this knowledge of w ~ ^ ~ 0(-\/ln 1/x), which is the change of the exponent 
of Ux, affect spreading, which, naively, one also would expect to follow this trend? 

The key lies in Eq. ( pQ), used outside the target. Take the noise term as an example: There the Ux in £xy 
and elsewhere in Eq. (pG) is 1. However, when y is inside the target while integrating, one gets a contribution 
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that changes Ux away from 1 via ( |60| ) since Uy is nontrivial. Imagine the target is characterized by some 
radius parameter, then this integral over y would be cut off there as 1 — [/^ is different from zero only inside: 

^1 

^^^fe^ (62) 

^0 if \y\>R 

Clearly, here the I — Uy factor determines the region of non vanishing random noise contributing to the 
evolution. It is a simple exercise to see that the force term & follows the same pattern. 

All these facts combine into a scaling argument: using the natural scaling properties of the ingredie nts in (|60| ) 
allows to scale the size of ^ ~ 0(-\/ln l/x) out of the equation almost everywhere using x — > y^ln 1/a; x, 
y -^In l/x Everywhere, except in the radius parameter which is now the only ingredient affected by 
the typical growth rate according to i? ^ ^ylnl/x R. This argument does not at all depend on the way such 
a radius parameter is introduced, it only relies on the fact that the only scale invariance breaking feature 
involved is the size of the target. Note, however, that it refers to the area of activity only -the area where 
Uy is nontrivial- not the measured cross section. These two notions become equal only when the blackness 
limit has been reached. 

The picture of evolution then is the following: During the whole evolution, the systems area of activity 
is growing like In l/x, the measurable cross section, however, will only start to follow this trend once the 
blackness limit has been reached. From then on In l/x type unitarization is visible in the cross sections, 
with little to distinguish a proton from uranium targets up to a global rcscaling with the cross sections at 
the "unitarization limit." One thus expects scaling laws that relate the evolution of different targets in this 
region. Unitarization is thus viewed as the limit in which the area over which the system appears to be black 
catches up with the total area, which scales with In l/x no matter which x one looks at. This is a remarkably 
concrete and hands on scenario to emerge from the rather mathematical elegance of Eqs. (H) and (||). 

The focus for further study then will become the question for saturation "efficiency," the speed at which an 
initially grey system becomes black on the one hand and the relative size of not yet saturated tails on the 
other. The negligibility of such tails is the main approximation underlying any scaling laws referred to above. 

The present equations reduce to a BFKL scenario at small densities and therefore imply power like saturation 
"efficiencies" when starting from moderate x. Note that the actual location of the unitarization limit will 
depend on target properties. As one is not talking about the asymptotic region of the BFKL equation -quite 
to the contrary- different components of the wave function will grow with different speed and it will depend 
on the amount of overlap with fast growing eigenfunctions if the unitarization limit is reached sooner or later. 

To turn back from general features to prospective numerical simulations, there is one item left to be supplied: 
an initial condition for the evolution, which will have to be adapted to the physical process to be explored. 

What one has to do is to generate an initial ensemble of configurations that, when averaged over, yields 
some physical cross section. At this point there are a few possibilities readily available. The first is to just 



follow 1 21 1 as closely as possible. This would limit the discussion to DIS per sc. The second somewhat more 
sophisticated approach would start with the McLerran-Venugopalan (MV) model for large nuclei H^,^. This 
also has something to tell about about the dipole cross section for fixed, small x as rudimentarily laid out 
in [po[ , but is more generic. This approach, however, should be updated to incorporate overall color neutrality 
as suggested in [Q. Although both approaches -the present RG and the MV model- are formulated with 
the help of statistical weights, the variables are not the same. The MV model is formulated in terms of what 
essentially corresponds to the exponent of U, the variables of this RG. This needs to be adapted. 



^^The scaling Him pnsj nn of ^ has to be length ^ , c.f. Eq. (|6l[). This is also consistent with having a scaling dimension for U 
of zero through Eq. (KOI) 
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8 Implications: correlators at vanishing x, universality and scaling 
laws exemplified 



The first and most obvious consequence of of the known form of Z at vanislring x, Eq. ( p4| ) inside the target, 
the area of activity, and Eq. outside, allows to make generic statements about the correlators in the 
"unitarization region," provided that the influence of tails that have not reached the saturation limit is small. 
This seems to be borne out by simulations in the large Nc limit One would also expect this to be of 

less importance with increasing target size as the system there starts from a situation with more gluons 
involved even initially and the tails proportionately smaller. 

One finds then that generic correlators with all distinct coordinates vanish inside the target and become 
equally trivial outside. Knowing that the area of activity will grow logarithmically as parametrized by a 
radius parameter one may write in the saturation region 

([/« L/W), ^0(?(any of{|a;i|, . . . , |a;„|} < Rj,{x)) 

+ 1 (g) . . . ® 1 0(an of {|a;i|, . . . , |a;„|} > 

iff all the Xi are distinct. Here the ""Q functions" are one if the condition indicated is met and zero otherwise. 
Exceptions to the above occur if one allows to have factors of U and C/^ at the same point as for example in 

+ e{\x\,\y\ > RA{x))5,jSki , 

due to the UxU]^ ~ 1 constraint in the measure. Similarly, the other constraint in the Haar measure, 
det J7 = 1, induces nontrivial expectation values with exact multiples of Nc factors of ?7s or [/^s inside. For 
example 

([C/ccJjiji ® ■■■® [C/x„Jj„^j„Jb -^Cx, eii...iiv,eji...iiv,^(a'i < ^.4(2;)) 

(65) 

+ ^im ■ • ■'5iN,iiv,^(all of{a;i, . . . ,a;jvj > RAi^A) ■ 

Any improvement beyond the naive step function scenario would naturally describe delayed saturation in the 
transverse tails of the wave function and would lead to replacing the step functions by some other effective 
"profile function" of the target. 

With this proviso one may calculate the asymptotics of specific correlators such as the dipole cross section in 
DIS; (tr(l — UxUy))b- The impact parameter integral will be cut off by the "outside" term in Eq. (^) and 
one finds 

d^ix + y) (tr(l - UxUy)^, ^ fA{TrR\{x))Nc{l - S^^h . (66) 



This is obviously universality of the dipole cross section, but not "saturation" , a feature also observed in the 
large Nc limit as reported in Ref. ||2^. The lack of saturation, of course, is a direct consequence of the global 
Inl/cc growth of the active area observed above. 

Here /^(vri?^) has the interpretation of the integral over a profile function for the target, the knowledge 
of which depends on whatever one is able to correctly implement when writing the initial condition for the 
evolution. To find the total cross section or the F2 structure function, this is to be convoluted with the 
correct photon wave function factors, but even without doing so, one more striking consequence is obvious: 

The above expression is only sensitive to the value of this wave function at its large effective momenta via 

('2') 

the 5xy and vanishing effective momenta via the 1 term. This leads to a complete loss of information about 
the momentum transfer of DIS, , in the limit of vanishing x. 

To come back to the scaling argument, one should expect that cross sections such as the above scale with 
their respective integrated profile functions /.A(7ri?^) (or equivalently with their ratios at the initial value 
of x). X evolution should be universal with this scaling law. Interestingly, this feature is already present in 
the large Nc limit, as demonstrated numerically in pl[ | for the dipole cross section in DIS. It it worth noting 
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that the simple scahng behavior found there according to the ratio of areas R\/ R% may be a consequence 
of the simple scaling behavior put into the initial conditions. More sophisticated approaches based on the 
McLerran-Venugopalan model [^E^, as used in ||2^ may lead to more complicated scaling factors with A 
dependence varying from A to A?!^. It is intriguing to note that in this numerical study scaling seems to 
hold everywhere, implying that saturation is very efficient and tails do not play an important role. It deserves 
further investigation to find out if this is due to the initial conditions chosen there or a rather more generic 
feature of these evolution equations. 

9 Conclusions 

The bridge spanned by this paper, from the presentation of an infinite hierarchy of evolution equations given 
originally in Ref. j^, summarized here in a single equation for a generating functional, to the functional 
Fokker-Planck or, more precisely, Brownian motion and Langevin representations has led one to a point from 
which the whole content of this "small x effective field theory" is much easier to understand. 

This has allowed to prove the existence and attractiveness of a fixed point at small x using simple arguments 
like positivity of the associated Fokker Planck Hamiltonian. This on its own makes it possible that there 
are universality features in experimental cross sections as one goes to smaller x. The nature of the evolution 
mechanism, including the driving force behind the unitarization mechanism have become apparent at this 
stage. 

In particular the Langevin equation has proven to be a formulation that encodes much of the physics intuition 
about what should be going on in this dynamical regime in formulae in a one to one correspondence. 

It has allowed to describe the unitarization mechanism at small x in very physical terms as the limit in 
which the area over which the system appears to be black catches up with the total area, which itself scales 
with Inl/a; throughout evolution. The mechanism driving the race towards this limit may depend on the 
observable at hand, for the DIS cross section it has been demonstrated that it is the BFKL evolution that 
drives the system towards unitarization, quite in agreement with standard expectations. 

With such a clear distinction between the onset of unitarization and the strict mathematical fixed point of 
the Brownian process itself, with the dynamical variables completely filling the configuration space, the latter 
has been found to be of only limited interest, as it refers to the limit at infinite energy with equally infinite 
total area and cross section. It is the much larger and experimentally accessible saturation region with its 
In 1 jx growth of the cross section that is the physics signal to be looked for. 

The results obtained for correlators in the saturation region clearly exemplify how the system loses information 
about the initial conditions already at this stage and give a clear idea about the kind of universality to be 
expected at small x. 

It is the yTnT/x erasing behavior of the Brownian process that drives and characterizes these aspects of 
the evolution. But there is also conceptual clarity encoded in this formulation: Any such stochastic process 
is accompanied by a loss of information, a fact nicely in correspondence with the Wilson RG interpretation 
underlying the integrating out of degrees of freedom used to (re)derive the initial version of the equations in 
this paper. 

At the same time the effective field theory nature encoded in the evolution equations has been complemented 
with a method that opens the doors to numerical simulations. The Langevin equations derived are suitable 
for this task and should help answer many if not all of the remaining "dynamical" questions. First on the list 
there is the question of saturation efficiencies and target dependence. This should give quantitative backing 
to expectations based on the BFKL limit and clarify open issues such as saturation in the transverse tails of 
the targets wave function as well as the amount of target dependence. 

There are of course farther ranging options. One of them is the possibility to use the results of such evolution 
as an input to describe the rapidity dependence of nuclear collisions in an extension of the work laid out 
in Q and made feasible numerically in [^,|2^. This would constitute one step deeper into the hierarchy of 
complexity detailed in Sec. |^ 
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Another intriguing question is that for the apphcation of the methods developed here to the JKLW equation. 
With the basic correspondence of the approaches known, an at least numerical Langevin type analysis should 
be possible. To do so, an analysis of the force terms, which need not be vanishing there, and possibly a 
reformulation in terms of eikonal operators are to be performed. This, however, needs more effort and space 
than available here and will have to be left for the future. 
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A Feynman rules 

Schematic Feynman rules used for these diagrams are given in the table below. All algebraic expressions 
needed in calculations are given in the main text. 



{5Ax5Ay)sA[b], the full gluon prop- 
agator in the presence of a "back- 
ground" field b 


Tr(r(r(r(^^b^innr(p 


Free gluon propagator 




a quark line or U factor along x~ 




an antiquark line or U'' factor along 
x~ 




The location of the color field (G"*"') 
of the target in x" ordered diagrams 
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B The vanishing of force terms 



In this appendix, the vanishing of the force terms in the FP equation is demonstrated by looking at the q 
component in somewhat greater detail. There one finds 



1 

2^ 



ICxzx (-2 Nc [Uxh + 2 [C/.]y tv{UiUx)) 
d^ylCx.yN, (-2 [Uxh + [UxUp,]^^ + [U.UlUx]^^ <5(2J 

ICx.x {[U.lj tr{UlUx) - Nc [Uxlj) \ 



(67) 



= 



A similar cancellation occurs in the q component. Note how the divergent terms cancel among the derivatives 
of X before the finite remainder cancels against a. 



C Transforming the evolution kernel 



C.l Lie derivatives are the natural objects to consider 



Suppose you never heard of Lie derivatives before. Do they come up naturally? They do: The first step is 
to express derivatives "along" SU{N, 
wants to respect the constraint UU^ 



to express derivatives "along" SU{Nc) via the independent and -g^. This is done observing that, if one 



1, one needs to replace for instance by 



SU 



_5_ 



kl 



which clearly does the trick 



r (Mc/t 
\Iwm 



- WMWU)af3 



uesu{N^) 



(68) 



(69) 



on SU{Nc). The only thing missing to turn the above into Lie derivatives is that they do not obey the correct 
commutation rules. This, however, can be achieved for the price of a variable change. Proper Lie derivatives 
on SU{Nc) may be written just like above, but with U'^'^'^ dependent AI. For this example letting M Ut^, 
one finds 



M, 



5_ 

'SUij 



5 



With A = a £ {1, . . . , A^^ — 1} one easily extracts the SU{Nc) commutation relations for [Ut'^]ij-^.. 



(70) 



^1 wn'ij 



for U £ SU (Nc) while still taking the derivatives with respect to U and W independently. 
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C.2 A little algebra 

There is a little catch that one has to deal with when trying to recxprcss things in these terms which has 
already been anticipated in the above notation. In order to be able to represent all of S/6uj^ps components 
one has to allow A to run over A = {0, . . . -1} and supply i° = l/^/57Vc 1, such that tr(t'^i^) = 
for all A, B € {0, — 1}- To perform the transformation, one writes 

and one is left with inverting ti-{t^Ut^) and ti{t^t^U'^). Using completeness 



one gets 

2tr(t^C/i'^) 2iY{t^U-h^) = 2tiit'^UU-h^) = S'^^ 
2iY{t'^t^U^) 2iY{t^f{U^)-^) = 5'^'^ . 

The transformation is then done according to cither variant in 

Mij = 2iY{t'^U-^M)[Ut%j ^ -2tr((C/t)-4'^A/)[-i^C/^]y , (71) 

taking proper account of any issues of chain rule when used for the Lie derivatives. Note that here "—1" is 
really "—1" as in "inverse" not "f ." 

To see how ( [n]) affects the FP equation, explore how x transforms. When doing so, chain rule terms appear 
in combinations like 

^M.,.^(2tr(t'^^-M)[^t1,) 

= [Ut%, — {2iY{t'^U-Hl))+NMU~^M) (72a) 



or 

-Mi. 



-{2tr{{U^)-h^M)[t^U^, 



= [-t''U%^{-2tv{{U^)-h^M)) + NMiU^y^M) (72b) 
where I have used -s^[Ut%j = 5^k5-jij^[Ut%i = iVctr(t'^) = 5^^Nc^/n7J2 to get the last term. 
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C.3 Results for the evolution kernel 

Chain rule aside for a moment, the components of the transformed x a-rc then of the form 

as2[t''U-%[x%]^^J^2[{Ul)-h%,) = -4^^ Jd'-z JC^^y (73a) 



X {tv{t''U^^U.t^UlUx) + tv{t^Ul{Ul)-h^UlU.) - tr(tV) - tv{t''U-\ul)-^t^UlUx) ^ 
a.{-2[{Ui)-h%,)[x%]^^^^2[t^Uy\u = -4^ jd^z Ky.^ (73b) 

X {ir{t^ulUyt^Uy'u.) + tv{{ui)-h^uiu.t^ul) - ir[tH^) - ivmrh^ulUyt^u;')) 

aalt^U^^ {xl%\^ u^\^^Vy\^ = -4^^ ^d^z K.^,y (73c) 
X (iTit'^U^^U^t^UlU^) + Irit'^UlUyt^Uy^U^) - tiit'^t^) - \.r{t''U^^Uyt^Uy^U^)^ 

aMUi)-'t%,[x%],^^,2[{Ul)-H%r, = -4^ jd^z IC^^y (73d) 
X (tr(t'^[/t([/t)-4'^[/tC/.) + trmrh^UlUJ^Ul) tr(([/t )-4'^^t (^t)-i^Bf^t) _ tr(t^i'^)' 
The key features of these expressions are as follows: 

• in SU{Nc) all of them vanish if either A or B is 0. Only octet components survive. All matrices 
involved now are explicitly in the adjoint representation. This does not hold outside SU{Nc) for the 
terms involving Us at z. 

• Combining Eqs. ^) with Eqs. ^) one concludes that all chain rule terms also vanish in SU{Nc). 

• in SU{Nc) all four of the above are equal. 

• This unique common form even factorizes naturally into the square of a much simpler factor. It will be 
denoted by xSy[C^] and reads 

asXlyiU] Jd'z IC^.y{u-'Uy + U^^U, - U^^Uy -iy\ (74) 

The factorized version is given in Eq. (pi]). 

D Reduction to the gauge group: details of proof 

It is actually instructive to see how Eq. (^) comes about. The cancellations involved will make it clear how 
evolution equations for correlators are recovered from the new form of the evolution equation. The easiest 
way to prove Eq. ( ^ ) is again by doing it for a generating functional for all those correlators in one go. In 
other words, one simply compares the action of the two operators on arbitrary correlators in the guise of 
e'^Bxtlt^.-'^l.pq The identity to prove then is simply 



Xuv c-ri- ^ 



2SUu SUy 

With the vanishing of the force term known, Eq. (^6|), the l.h.s. is straightforwardly given by 
1 <5 _ S 



(75) 

c/c/t=i 



A 



2SUu SUy 



(7C/+ 



(76) 
c/c/t=i 



^^Again, the arbitrariness of the correlators results from that of the sources J and the equations for correlators follow by 
differentiation w.r.t. at J = 0. 



26 



while the r.h.s is 
1 



_ 1 r 

(7(7+ = l 2 



4{ [OV"^^xf^)tr((jt);^C/^t'') +x^^^tr((jt)^C/^t-<^)5(2j 



(7C/t=l 



(7C/t = l 



=2 JT 



- iV 



(7(7+=l 



;7;7t=i 



q.e.d. 

Note how the ct terms emerge. Quite a surprising combination. 



(77) 



E From Lie derivatives to zero modes of the FP Hamiltonian 

E.l Lie derivatives and transverse local gauge transformations 

There is more to be said besides^ 



(78) 



to characterize the Lie derivatives as an infinitesimal gauge transformation. They truly do implement trans- 
verse local gauge transformations, a feature that can be summarized by 



c^S^v^Ja,7^(L/))^ = 7^(L/,X7^(r)(afCj) 
c.s^v^Ja.7^(c/-l))^ = -^S7^(^'^)7^(c/-l)(9f4'J) • 



(79) 



This implements an infinitesimal change to U in such a way that (restricting to the fundamental representation 
for brevity) 



(80) 



Clearly this is an infinitesimal local gauge transformation (with respect to the transverse coordinate) of a 
pure gauge potential hi = {diU) ^. 

This is in one to one correspondence with the Langevin equation which also implements an infinitesimal 
transverse local transformation according to 



in a (randomly) Inl/x dependent manner in each step of the evolution. 

In summary, all the evolution does in a finite x interval is to evolve any Ux{yo) into 



(81) 



(82) 



'^'R,{Uy) indicates Uy in a given representation TZ. Just leave it off for tlie fundamental representation. 
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where y = h\\ / x is the rapidity. 

It is important to note, however, that this notion of a local gauge transformation only applies to the transverse 
coordinate dependence. Any y dependence is related to coarse graining in the ± plane and will not behave 
like a gauge transformation with respect to these coordinates. 



E.2 A unique zero mode 



With the nature of the evolution clearly established, the zero modes of the FP Hamiltonian arc also quite easily 
listed: Any function of only scalar, (transversally) locally gauge invariant operators will do as a candidate 
for Zoo[U]. The list of such arguments is in practice much shorter that it might initially appear. Indeed, 
purely (transversally) 2-d operators are generically trivial: There the only ones one could write would have 
to contain at least one factor of Fij\U~^dU] and that is, of course, zero. 

All other apparent candidates would necessarily involve a ± derivative, say via i^"*"', and will necessarily not 
be invariant .p*^ In particular the weight of the MV model is not an invariant of this RG, neither in its initial 
form nor its "color neural" setup according to |^^. The same argument excludes any "natural" topological 
term of the Wcss Zumino type / e^^^tx[U~^[dAU)U~^[dBU)U~^[dcU))vi\ which the third coordinate would 
have to be taken again from within the ± plane. This exhausts the list of all possible candidates of any 
dimension and leaves a "trivial" but unique solution: Z^o \U] ~ 1 . 
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